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ABSTRACT 


The operational effectiveness of airborne laser hydrography 
systems, considering the optical environment of the coastal waters 
of Oregon, Washington, and the Gulf Coast states, is examined. 

The best times of the year are predicted for conducting laser 
meenymetry, considering the temporal and spatial variability of 
u properties due to seasonal effects, and charts of season” 
ally averaged optical measurements are given. Original formulas 
to convert beam attenuation coefficients and Secchi depth measure- 
ments to irradiance attenuation coefficients are included. The 
number of irradiance attenuation lengths to the bottom depth (Kd) 
are used as the indicator to estimate areas where laser hydrog- 
raphy systems would be successful and are shown by season and 
region. The conclusionsof this thesis are that airborne laser 
meerography is not practical in the coastal waters of Oregon 

and Washington, would be practical in limited areas of the western 
Gulf Coast, and would be very practical in the eastern Gulf Coast 


area. Along the eastern Gulf Coast a 38,800 nmi~ area, delineated 


by a Kd = 4 contour, is judged surveyable by laser. 
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I. INTRODUCTION 


A. HYDROGRAPHY AND THE DEVELOPMENT OF LASER BATHYMETRY 

Hydrography has been defined as, "that science which deals 
with the measurement and description of the physical features 
of the oceans, seas, lakes, rivers, and other waters, and their 
adjoining coastal areas, with particular reference to their use 
for navigational purposes" (U.S. Naval Oceanographic Office, 
1966). Hydrographic surveys are used to produce nautical charts 
and related information which satisfy navigational, engineering 
and marine scientific needs and contribute to national goals 
such as ocean resource management and national defense. However, 
these products may become quickly outdated, and new surveys are 
required because of changes caused by such natural processes as 
winds, tides, earthquakes and because of man-made changes result- 
ing from construction of ports, channels, breakwaters, and 
Pipelines. New surveys might also be needed because of changing 
requirements, for example those associated with deep draft 
tankers. Since the 1930s most hydrographic surveys have been 
performed by ships equipped with sonic depth sounders to obtain 
bathymetric data. 

A new method for hydrographic surveying, airborne laser 
bathymetry, was shown to be a fast reliable technique to obtain 
Macnymetric data by the U.S. Naval Oceanographic Office (NAVOCEANO) 
Been their Pulsed Light Airborne Depth Sounder (PLADS) in 1969 
(Bright, 1973). The PLADS system used a frequency-doubled 


Neodymium Yttrium Aluminum Garnet (Nd:YAG) laser (Rattman and 
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Cunningham, 1969). NAVCCEANO's interest in such a system 
resulted from immediate requirements stemming from the Vietnam 
War for better charts of shallow coastal waters within range 
Brzenemy controlled territory (Bright, 1973). 

Since that time further work has been done by NAVOCEANO, 
NASA, National Ocean Survey (NOS), and the Australian Weapons 
Research Establishment. In 1974 NAVOCEANO tested a NASA-owned 
neon laser and in 1975 obtained a frequency-doubled Nd:YAG 
E tor their Coastal Aerial Photo-Laser Survey (CAPS) system 
(Crandall, 1976). In 1977 NASA,with DOD and NOS sponsorship, 
Ma testing their Airborne Oceanographic Lidar (AOL) as a 
hydrographic data acquisition system (Guenther and Enabnit, 
EN Ihe success of these tests influenced the Naval 
Beeanographic Research and Development Agency (NORDA) to place 
meer Hydrographic Airborne Laser Sounder (HALS) in the procure- 
ment stage in 1978 and influenced NOS to implement a development 
plan for a more sophisticated system in 1979. The Australian 
Weapons Research Establishment has built two systems, one for 
research and one for operational surveys (National Ocean Survey, 


m9). 


B. MAJOR BENEFITS OBTAINABLE FROM LASER BATHYMETRY 

The National Ocean Survey (1979) has evaluated airborne 
laser hydrography and has found four major benefits realizable 
over ship/sonar hydrography with this technique: cost savings, 
manpower savings, capability of increased production, and im- 


provement in the quality of marine charts. 
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The NOS (1979) study showed that projected cost savings 
would be achieved through the speed with which an area could 
be surveyed. Figure 1 shows the costs per unit area as a 
function of area surveyed annually. The three curves represent 
amortized nonrecurring capital cost, operating cost, and total 
cost. The constant line approximates the cost of launch/sonar 
hydrography at $2730 per square nautical mile. The comparison 
of the projected laser cost per unit area of $438/nni* to the 
launch/sonar cost per unit area of $2730/nmi* For Есу 
utilized system, 2000 nmi^/yr, indicates that laser surveys 
would cost one-sixth of sonar surveys. 

The NOS (1979) study also examined manpower and production. 
NOS concluded that manpower effectiveness of laser over launch/ 
sonar would be five-to-one. Projected production for each NOS 
airborne laser system would be approximately 2000 пті” рег 
annum. 

Finally, the NOS (1969) study suggested that the quality 
of marine charts would be improved by laser systems because of 
Significantly increased spatial density and greater uniformity 
of distribution of soundings compared to sonar. The proposed 
NOS laser system was to collect 400-600 depth measurements per 
second with an average distribution of one per 20 mé, yielding 
measurements 4.5 m apart in all directions. This would be 
300 times the number of soundings per unit area of typical 
ship/sonar surveys. The increased density and more uniform 


distribution would provide a more representative chart. 
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РЕПЕРОСЕ ОҒ THE THESIS 

NAVOCEANO and NASA have shown that laser bathymetry systems 
may be practical and reliable. NOS has shown that major bene- 
fits for hydrography may be possible with a laser system. The 
purpose of this thesis was to determine the most effective use 
of airborne laser hydrography systems considering the marine 
environment. Specifically, the questions to be answered are: 
which areas of U.S. coastal waters may be most effectively 
surveyed using airborne laser bathymetry and at what times of 
the year? This study was confined to the Gulf and West Coasts 
because major studies have already been conducted by NOS for 


the East Coast (Enabnit, 1979). 
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ТІ. ATTENUATION OF LASER BEAM POWER AS THE PARAMETER GOVERNING 
ШЕН EFFECTIVENESS OF LASER BATHYMETRY 

A. INTRODUCTION 

Laser penetration of ocean water is dependent on the optical 
properties of sea water. This thesis examines the optical prop- 
erties which can be used to delineate those ocean areas which 
might be surveyed advantageously by laser bathymetry. Definitions 
and terms used to describe optical properties of sea water are 
those recommended by the Committee on Radiant Energy in the Sea 


given by Jerlov (1976). 


B. OPTICAL CLASSIFICATION OF OCEAN WATERS 

Jerlov (1976) has classified ocean waters in terms of the 
spectral transmittance of downward irradiance at high solar 
altitudes (Fig. 2). The Jerlov coastal water types 1-9 are 
characterized by increasingly higher amounts of yellow substance. 
Water types of decreasing irradiance transmittance indicate a 
spectral shift in the transmittance maximum toward longer wave- 
lengths. Selective absorption by particles and yellow substance 
causes greater absorption at the shorter wavelengths and shifts 
the transmittance maximum from 470 nm, the blue region, for 
clear ocean type I water to 550 nm, the green region, for 
coastal type 7 water. Figure 3 is a similar graph of Jerlov's 
water types but gives instead of transmittance the irradiance 
attenuation coefficient (K) for downwelling daylight as a 


function of wavelength. 
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С. TRANSMISSION WINDOWS AND LASER TYPES 

The preference for a particular type of laser to be used 
for laser bathymetry depends upon the wavelength of maximum 
transmission. The transmission window for Jerlov's coastal 
water types 1-7 is approximately 510-580 nm (Figs. 2 and 3). 
Laser types that have operating wavelengths within that spectral 
band are: the frequency-doubled Nd:YAG laser (used in the 
proposed NORDA HALS and NOS systems), which operates at 532 nm; 
the neon laser (used in NASA's AOL), which operates at 540 nm; 
the argon laser, which operates at 514 nm; and the dye laser, 
which is tunable over the blue-green spectrum (Ferguson, 1975). 
The frequency-doubled Nd:YAG is preferred by NORDA and NOS for 
two primary reasons: (1) its high peak pulse power and high 


pulse rates, and (2) its small size and weight (NORDA, 1978). 


ши НЕ PRINCIPLE OF LASER BATHYMETRY 

The technique of using a pulsed laser to measure water 
depths remotely can be explained with the aid of Figure 4. A 
short pulse of light is emitted from an airborne laser. The 
pulse of energy travels at the velocity of light and impinges 
upon the surface of the water. Approximately 3$ of this energy 
is specularly reflected from the air/water interface and inter- 
cepted by the aircraft receiver. Half of the time difference 
|between the initial laser pulse and this surface reflection 
yields the aircraft altitude. The remaining 97$ of the laser 
energy is transmitted into the water, wherein its velocity is 
decreased by about 25$. (The transmission across the air/water 
interface is dependent on the angle of incidence, on polariza- 
tion, and on sea state conditions. The value of 97$ is an 
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average value for angles of incidence between 0° and 45° in 
calm seas [Witt, 19791.) In addition the signal is 
exponentially attenuated by absorption and scattering within 
the water column. If the signal is of sufficient intensity, 
it will be reflected from the water/sediment interface and be 
detected at the receiver. The time difference between the 
surface and bottom-sediment reflections is used to determine 
the water depth. 
EV "TE: IMPORTANCE OF THE SYSTEM ATTENUATION COEFFICIENT IN THE 
LASER SIGNAL EQUATION 
A signal equation for a laser pulse transmitted from an 
airborne platform to the ocean floor may be given as follows 


(Avco Everett Research Laboratory, Inc., 1975): 


-2a һ-2үа 





: Py R(l-p)e 1 А тоң 

Pra = 
т.) п? 
n 
Bea = Received peak power [W] 
p. = Transmitted peak power [W] 
R = Bottom reflectivity 
О = Surface reflectivity 
a] = Atmospheric attenuation [km 5] 
= System attenuation coefficient of sea water (m +] 

А = Area of collector [m°] 
NS - System efficiency 
h = Aircraft altitude [km] 
d - Water depth [πι] 
n = Index of refraction of sea water 
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This equation shows that the received signal is a function of 
environmental and system parameters. However, for a given appro- 
priately designed and operated laser system, the major operational 
limitation is due to the environmental parameters. A simplified 


signal equation that indicates this is: 


P SR η ¿O 


EB 


2 


) / [T (h+a/n) 2 n 


where S = (Р А 1 


c SB 
For an aircraft altitude of 609 m (2000 ft), a water depth 

of 10 m, and representative values for the system and environ- 

mental parameters PL = 30 kW, R = 15%, p = 25, a, = O m 

y = 0.2/m, Ας = 0.073 eN Top = 1٦۔278‎ rana n = 1.33 (Avco Everett 

Research Laboratory, Inc., 1975) the following values are 

obtained: 


τ. UR WR oS 1-0 = 0.98, е 2918 - 0.86, 


Ор Е τα κ 10 W. 


The system parameter, S, 15 fixed for a given system or survey 
operation. Therefore, the ultimate environmental parameter 15 
the product of the system attenuation coefficient and water 


depth (yd) since it dominates all other terns. 


F. THE IRRADIANCE ATTENUTATION COEFFICIENT AS AN APPROXIMATION 
Witt (1979) showed that the system attenuation co- 
efficient (y) may best be approximated by the irradiance attenua- 
tion coefficient (K), for downwelling light, an apparent optical 
oceanographic property readily measured in situ. This approxima- 


MON was further confirmed by Krumboltz (1979) in a series of 
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tests for NORDA's HALS system. The coefficient K best approxi- 
mates y because after traversing several attenuation lengths, 
the photons of a laser beam undergo multiple scattering, 
spreading the shape of the beam and ultimately giving it an 
asymptotic radiance distribution equivalent in shape to the 
irradiance distribution of downwelling daylight. 

TS (1979) and Krumboltz (1979) have shown that 
Ever field of view (FOV) and altitude affect the accuracy 
of this approximation. Figure 5 shows the variation of y as 
a function of FOV and receiver altitude at one test site. 

At a receiver altitude of 500 ft, y = K with a FOV of 80 mrad, 
et 1500 ft, Y ~ K with a FOV of 40 mrad. The better approxi- 
mation at higher receiver altitudes is due to the greater laser 
divergence with distance. The NASA AOL system had a FOV of 
59-20 mrad (Avco Everett Research Laboratory, Inc., 1975), the 
NORDA HALS system a FOV of 0-30 mrad (Naval Oceanographic 
Research and Development Agency, 1978), and the proposed NOS 
system a FOV of 0-50 mrad (Avco Everett Research Laboratory, 
Inc., 1978). Figure 5 indicates that for the NASA AOL, the 
NORDA HALS, and the projected NOS systems the approximation 

1s good for airborne laser bathymetry systems operating at 

500 ft or higher but should be used with discretion at lower 
altitudes. 

G. IRRADIANCE ATTENUATION LENGTHS AND PERFORMANCE OF LASER 

BerHYMETRY SYSTEMS 

The performance of laser bathymetry systems may be specified 
using irradiance attenuation lengths, L, which are the reciprocals 


 Екайтіапсе attenuation coefficients, i.e., L = 1/K, the 
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distance in which the irradiance decreases by a factor l/e. 
The number of irradiance attenuation lengths to the bottom 
depth (а/т, = Ка) that various laser bathymetry systems can 
attain or are projected to attain is shown in Fig. 6. 

The NASA AOL system has achieved two attenuation lengths 
in the daytime and three at night (Enabnit, 1979), the NORDA 
HALS system was contracted to achieve 3.2 (Houck, 1979), and 
the proposed NOS system has a goal of four attenuation lengths 
(Enabnit, 1979). A ΠΠ. value can be attained at night due 
to the lack of background noise caused by sunlight. It must 
be noted that in most studies of laser bathymetry, beam attenua- 
tion lengths have been used in describing the capability of 
laser systems, even though Krumboltz 
(1979) has shown that the system attenuation coefficient was 
better approximated by the irradiance attenuation coefficient. 
The problem resulting from the use of beam attenuation lengths 
is apparent from Figure 5.) Thus some investigators have 
reported attenuation lengths five times greater than those given 


here (Guenther and Enabnit, 1978). 


H. NATURE OF THE PROBLEM 
The first thesis objective was to determine where in U.S. 
coastal waters laser bathymetry can be effective, using the 
irradiance attenuation length as the best indicator of the 
expected performance of airborne laser bathymetry systems. 
The second objective was to predict the best times of the 
year for using laser bathymetry. The optical properties of 


coastal waters may have temporal and spatial variability due 
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to seasonal changes in currents, winds, runoff, upwelling and 
other physical processes, as well as to biological activity. 
Other effects of shorter duration and more limited extent are 
storms, plankton blooms, tides, and man-made pollutants. Only 
the temporal and spatial variability due to seasonal effects 

were investigated. These are summarized by figures of seasonally 


Slotted optical data. 
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III. DATA COLLECTION AND PROCESSING 


A. LIMITATIONS OF MARINE OPTICAL DATA IN COASTAL WATERS 

One of the major tasks in this study was the collection of 
Marine optical data. Optical data for coastal waters are 
scarce; this is especially true for irradiance attenuation data. 
The most common type of optical measurement, the Secchi depth 
(2), is usually obtained by marine biologists studying the 
euphotic zone. The next most common type of optical measure- 
ment for the Gulf Coast area is beam attenuation.  Irradiance 
attenuation measurements, the optical data actually sought, 
are almost nonexistent for U.S. coastal waters. This scarcity 
of irradiance attenuation data required the use of Secchi depths 
and beam attenuation data. Other marine optical measurements 
such as volume scattering and turbidity measurements such as 
those recorded in Jackson or Formazin turbidity units were not 
considered. The in situ data sources used are listed in 
Appendix A. A limiting condition on all optical data used was 
that the water depth at the ocean station where the optical 
measurements were made had to be less than 200 meters. Since 
laser bathymetry systems will probably be limited to approxi- 
S IV four irradiance attenuation lengths (Enabnit, 1979), 


0 п depth for a K of 0.05 m! 


), 200 meters was considered 
Mere than sufficient. 

Satellite data may be indicative of the optical quality of 
coastal waters. Therefore, selected Landsat and NOAA-3 


satellite imagery were qualitatively examined for turbidity 


Patterns that correlated with in situ beam attenuation data. 
22 





Another necessary task was to find conversion formulas to 
convert Secchi depths (2_) and beam attenuation coefficients (C) 
to irradiance attenuation coefficients (K), because irradiance 
attenuation length was the chosen operational parameter. A 
related problem arose due to the spectral dependency of optical 
measurements. The optical data would be most useful if it were 
obtained for a wavelength of approximately 532 nm (wave- 

ШЕГЕН ОГ Тһе frequency doubled Nd:YAG laser), but such data 
are scarce. This scarcity required the use of optical data 
measured with broad bandwidths or at wavelengths other than 
532 nm and the formulation of the necessary conversion relations. 

The third major task was the establishment of a data base 
with common optical and measurement parameters. Data from sites 
within a certain ocean region and collected during the same 


optical season were averaged. 


B. DATA REDUCTION AND CONVERSION 
l. Irradiance Attenuation Coefficient (К) Data 
The irradiance attenuation coefficient, K(A), is an 
apparent, spectrally dependent optical property that measures 
the extent to which diffuse downwelling davlight diminishes 


exponentially with depth in water. Trus, 


H =нде F LA) 2 
z O 


solar irradiance at sea surfaces 


where: H 
O 


H 


> downwelling irradiance at depth 2 


Z = depth of measurement 


K(1) = irradiance attenuation coefficient* 





* 
Normally (A) will not be written, as it will generally be 
assumed that K is for a specific wavelength. 


“A 





Depth averaged values of irradiance attenuation co- 
efficients (K) were required. Data were usually obtained in 
percent transmissions ($T - 100 H,/H,) for a series of depths. 
TO obtain a depth-averaged value, K, for the entire water column 
the following equation was used in a linear regression using 


the method of least squares: 


ln(T,) = In(H,/H ), = -К2„ 


n = index of data palrs 
K in units of пі was the slope of the straight-line model with 
depths, 24! in meters as the independent variable and the corre- 
sponding transmissions, Т, as the dependent variable. 

Ideally, irradiance attenuation data used for this study 
should have been measured at A = 532 nm. However, few such data 
were found and irradiance data at other wavelengths had to be 
used. Most of these data were taken with irradiance meters 
having a photopic response with a peak transmission near 555 nm. 
The broad bandwidth of a photopic response coupled with the 
water's selective absorption by wavelength with increasing depth 
yields an effective response which is shifted towards the 
wavelength of maximum transmission for the given water. For 
this reason photopic K measurements were considered useful. 

Other K measurements used in this study were obtained 
with a quantum irradiance meter. A quantum meter is normally 
sensitive in the 350-700 nm bandwidth interval and has a spectral 
energy sensitivity directly proportional to wavelength (Jerlov, 
1976). To convert from quantum values, Ko to K the curves of 


Figs. 7 and 8 (Jerlov, 1976) were used. Given quantum transmission 
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measurements at a station, Fig. 7 was used to determine a water 
type. Figure 8 was then used to determine the к/к ratio: 
This ratio was then multiplied by the depth-averaged quantum 


irradiance attenuation coefficient (K..) to obtain the irradiance 


Q 
attenuation coefficient (K). Again the broad bandwidth of the 
guantum meter becomes narrowed to a green band withdepth in 
Epastal waters. 
2. Beam Attenuation Coefficient (C) Data 

The beam attenuation coefficient, C(A), is an inherent, 
spectrally dependent, optical property that characterizes the 
attenuation due to absorption and scattering by a collimated 
beam of monochromatic light traversing a fixed path-length of 
homogeneous water. Thus, 
e C Or 


FL = Ir 


where: Po initial radiant power from projector 


Е, = residual radiant power measured by receiver 
r = path length of measurement 
C(A) = beam attenuation coefficient* 


Depth averaged values of beam attenuation coefficients (С) 
were computed for this study. Data were usually expressed in 
percent transmission per meter, 3T = 100 FL/Foı for a series 
of depths. To obtain a depth-averaged C for a water column, 
merical integration of the transmittance profile with depth 
was performed. First, the water column at each observation 
station was divided into layers so that each beam attenuation 


measurement was at the center of that layer. Two exceptions to 


л. 


* 
Normally (X) will notbe written as it will generally be 
assumed that C is for a specific wavelength. 
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this procedure were the surface and bottom layers, because layer 
boundaries between measurements were kept equidistant. Second, 
each layer was assumed to contain homogeneous water, and the 
beam attenuation coefficient was calculated for each layer using 
the equation: 

1] 


С. =-1n(F,/F,) [m 


Third, the transmittance of the entire water column was calcu- 


lated using an iterative procedure: 


T, 7 T4, exp (-C, rj) 


where Cu = coefficient for the layer (m 1) in question 


r = width of the layer (m) in question 
T ss incoming transmission to the layer in question 


and outgoing transmission of the previous 


layer. 
т outgoing transmission of the layer in question 
and incoming transmission to the next layer. 
к= тоо» 


n = number of layers 


Fourth, the depth averaged C was calculated by using the equation: 


E 


“l/r, In (T7100) 


Es 


where С - depth averaged beam attenuation coefficient (m 


r water depth at the station (m) 


2 


T final Т, of previous equation 
Not all the beam attenuation data used for this study 
Were measured at y = 532 nm. Beam attenuation data obtained for 


Oregon coastal waters (Pak, 1979) were measured at 660 nm and 
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had to be converted. A procedure to convert this data to 532 nm 
was formulated using data collected by Zaneveld, et al. (1978) 
 ШМопгегеу Bay, California. They used a spectral beam trans- 
missometer that simultaneously measured at six wavelengths at 
various depths. From these beam attenuation data those for 
45 stations at three wavelengths (500, 550, and 650 nm) 
measured at six different depths (0-1, 5, 10, 20, 40, and 60 
meters) were selected for this study. Linear regression models 
were computed (x, y: y = mx + b) for the following sets of 
КОО ед: 1|С(650), С(500) |, |С(500), C(650)], [C(650), C(550)], 
and [C(550), C(650)] at each depth. Table I summarizes the 
procedure followed in computing a mean slope (m) and a mean y- 
intercept (b) to obtain the two equations: 

e 500) Fr 122€ (6509070442 


C (550) 


mao C (650) - 0.35 

The final equation, obtained by interpolating between the above 

two equations, was used to convert C(660) to C(532). 
۰۰ 6 7 

ге С(650): 0.425 < C(650) < 2.00 [m ] 

Beam attenuation coefficients were converted to irradiance 
attenuation coefficients using a relation computed by Shannon 
ERU) for C(532). 

9۹۹ ۹ ۹۶۹۶١۶ ؤ٦‎ 

۳۰ >۳ ںی ۷۷۷۹ء‎ 86 [m t] 


3。 Secchi Depth (22) Баға 


Most of the optical oceanographic measurements used were 
Secchi depths. This is the depth at which a white Secchi disc, 
usually 30 cm in diameter, lowered into the sea and viewed from 
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directly above, disappears from view. Unfortunately, the proce- 
dure has never really been standardized (Tyler, 1968; Williams, 
1968; Holmes, 1970). Factors affecting the visibility of a 
Secchi disc include solar altitude, cloud cover, sea surface 
reflection and refraction, ship shadow, water color, observer's 
visual acuity, and height above the water surface (Tyler, 1968). 
Thus, Secchi disc measurements are somewhat subjective and 
imprecise. 

As are the beam and irradiance attenuation measurements, 
the Secchi disc measurement is spectrally dependent. For the 
purpose of this study, the observers should have worn eye 
glasses with, for example, Wratten 61 filters (Williams, 1968). 
Again, such Spectrally homogeneous data are not available. 
Williams (1968) investigated the relative response of the 
photopic eye (555 nm peak response) with and without a Wratten 
#61 filter (dominant wavelength 530 nm). Because of the 
selective absorption with depth of ocean water, the errors 
generated in uSing the photopic Secchi measurements were not 
hE to be of significance for this study. 

Secchi depth measurements were converted to irradiance 
attenuation coefficients by methods based on geographic 
region. Formulas published by Poole and Atkins (1929), Graham 
(1966), Otto (1966), and Holmes (1970) were not considered ade- 
quate for this study because of poor statistical procedures or 
because the data were obtained in a different type of water. 
Instead, data from various sources were used to develop regional 


Or mulas to convert Z to K. 
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The conversion formula used for Eastern Pacific coastal 
waters was based on photopic irradiance and Secchi depth data 
published by Callaway and McGary (1959) and Holmes (1970). 

The data selected from Callaway and McGary were from 19 stations 
approximately 150 nautical miles off the U.S. coast and north of 
35°N latitude (Table II). Added to these were 13 stations in 
Goleta Bay near Santa Barbara, California, observed by Holmes. 
These 32 stations were used and processed in the following 
manner. First, the Callaway and McGary irradiance attenuation 
data were transformed to a depth-averaged value for the same 
water layer as the Secchi depth measurement. This value (К) 
was computed by using only those irradiance transmission values 
which were measured or interpolated at depths equal to or 
shallower than the Secchi depth. The procedure used to calcu- 
late К. from these transmission values was the same linear 
regression technique mentioned in Section III.B.l. The Holmes 
irradiance data were already in the correct form. Second, two 
linear regression models were used to compute conversion 
formulas. The first linear regression model used the Secchi 
depth reciprocal ο... as the independent variable and κ. as 


the dependent variable. The equation obtained was: 


К Тот АЩ ROUGE 


-1 
Е ) 


95% confidence interval 


m = 1.21 + 0.10 


O 
N 


0.06 + 0.02 


r = .98 
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In the second linear regression model δε was interchanged 


with κ... The equation obtained was: 


=l 


Z ο 76 K 7 9 7 


] 


95% confidence lnterval 


S 


m = 0.78 + 0.07 


b 


0504 55:50:02 

22 20.95 
The final conversion formula used is the mean of the above two 
equations: 


zi 


+0 257 + 5 [m 
S S 


where Z.: 75٣-٣ ٣٠32 тп) 


Figure 9 is a graph of the data points and the regression lines. 
The conversion formula used for the western Gulf Coast 
waters (Longitude > 89.51) was based on data published by 
Kamykowski, et al. (1978), quantum irradiance and Secchi depth 
data for one station off Texas observed on four different 
cruises. Eleven of the twelve published measurements were 
used. First, the data were transformed from quantum irradiance 
transmission to depth-averaged quantum irradiance attentuation 
coefficients for the same layer of water as the Secchi depth 
Measurements. Second, these coefficients were converted from 
the Ko to K by using Figs. 7 and 8 from Jerlov (1976). Table 
III summarizes the results. Third, linear regression models 


were computed resulting in the following equations: 


- _ -1 -1 
K_ = 1.142” +0.04 [m 


] 


95% confidence interval 


m = 1.14 + 0.39 
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р = 0.04 + 0.05 


I= ۱۱ 
-1 
| 


95% confidence interval 


2 = 0.73K. - 0.01 [m 


ω ου 0-25 
b = -0.01 + 0.05 
C =el 
The final conversion formula used is the mean of the above two 
equations: 
1 -1 


K. = 1.26 Z, + 0.03 [m 


] 


where 2: зе be S 26 [m] 


Figure 10 is a graph of the data points and the regression lines. 
The conversion formula used for the Eastern Gulf coastal 
waters (Longitude < 89.5°W) was computed by Shannon (Wite BE 51... 
1976) from optical measurements taken in the West Pacific and 
Eastern U.S. coastal waters. 
l 


ο 0.03 Та 
5 5 


-1 
| 
r= .95 


where Z.: 4 < 2. < 43 [m] 


The K data used were measured with two types of irradiance 

meters--photopic response or a peak response at 533.5 nm 

(Witt, 12979). This formula was determined by a linear re- 
-1 


gression model with 2. ` the independent variable and κ. the 


dependent variable (Shannon, 1979). The weakness of this 


zn 


equation was that 2. as the independent variable was assumed 


to be error free. 
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The Poole and Atkins (1929) data were examined due to 
the wide use of their formula, 
1 -1 
| 


Poole and Atkins used an irradiance meter with a photopic response 


K= l.7 Ze [m 


and made fourteen series of observations mostly at one station in 


the English Channel over the course of one year. An average K was 
calculated for the surface water layer (0-20 m) for each series, 


τ) --; 


Poole and Atkins then found the mean value of the fourteen < Z 


products to be 1.7. “he same Poole and Atkins data were used 


and this was multiplied by Z., giving the product K 


to recalculate a formula using linear regression and to verify 
the conversions of this thesis. The two computed linear regression 


equations are: 


p = 1.11 ο + 0.04 m +] 
95% confidence interval 
700 1 9 9 ο» 
ТЕЕ 
г = .90 
а оок 11 3٦پ‎ 


95% confidence interval 


m= 0.75 ЗЕ 0222 
b = -0.01 ο. ο ο 
r = 90 


The mean equation from the above two is: 
= - = -1 
Ko-20 = ص5‎ 2 + 0.03 [m | 
where 2: Bes Z. < 22 [m] 


Figure 11 is a graph of the data points and the regression lines. 
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The close agreement between these equations verified 
the conversion formulas used to estimate the irradiance 
attenuation coefficient from the Secchi depth. Table IV 
zes the regression equations, and Fig. 12 shows the 


mal conversion formulas. 


@ REMOTELY SENSED DATA 

Three types of satellites are able to provide data that 
indicate the optical quality of ocean waters: Landsat series, 
NOAA series and Nimbus 7. Selected imagery from Landsat and 
NOAA series were examined. 

Landsat data have been used with some success by numerous 
investigators studying currents tagged by suspended sediments, 
plumes, and dispersal patterns of suspended sediments (Jarman, 
MS; Pirie and Steller, 1973; Hunter, 1973; Erb, 1974; Carlson, 
8:157 1975; Maul and Gordon, 1975; Rouse and Coleman, 1976). 
The most successful studies have used image enhancement processes 
with in situ data simultaneously collected for calibration 
purposes. However, such enhancement was not available for 
this study, and only B&W imagery as received from the EROS data 
center of the U.S. Geological Survey were used. 

For the coastal region off Oregon MSS4 and MSS5 Landsat 
imagery taken during the March 11, March 28, May 12, and May 13, 
1975, overpasses were compared to beam attenuation data (1 = 
660 nm) collected on April 23-May 1, 1975 by Pak and Zaneveld 
(1977). No noticeable correlations were observed. 

For the region off Texas MSS5 Landsat imagery taken during 


the November 12-13, 1975, and May 28-29, 1976, overpasses were 
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compared to beam attenuation data (X = 528 nm) published by 
Berryhill, et al. (1976) and collected on November 15-21, 1975 
and May 21-25, 1976, respectively. Turbidity patterns were 
much more noticeable for the Texas imagery and some trends 
between the satellite imagery and the beam attenuation data 
were observed. These turbidity patterns observed in the 
andsat imagery were used as an aid in contouring the ات‎ 

beam attenuation coefficient data. Shideler (1979) extensively 
examined the data of all six cruises published in Berryhill, 

et al. (1976) and Berryhill, et al.(1978) with the related BsW 
Landsat imagery and indicated that turbid water masses through- 
out the inner shelf depths (less than 45 m) were apparent in 
the imagery. 

The NOAA-3 imagery for selected days throughout 1975 were 
used to observe upwelling off the U.S. West Coast. The relation- 
ship Of upwelling to seasonal variation of optical properties 
off the Northwest Coast is discussed in the next chapter. 

Meastal Zone Color Scanner (CZCS) data from the Nimbus 7 
Satellite were not available at the time of the writing of this 
thesis. However, it is mentioned here because of its enormous 
potential as an optical data source for the marine environment. 
Data can be obtained on magnetic tape, CRCST (Calibrated 
Radiance, Pigment, Diffuse Attenuation Coefficient, and Tempera- 
ture Tape), or film format. The irradiance (diffuse) attenuation 
coefficient is a computer parameter obtained on tape or film 
image. Since this scanner views a swath 1566 km wide, K data 
for enormous ocean areas could be obtained (Goddard Space Flight 


Center, 1978). 





Bee DATA PROCESSING 

Each measurement datum was keypunched onto a standard 80- 
character computer card along with its latitude, longitude, 
date, bottom depth (when available), measurement unit (meters, 
feet, etc.), converted K value, and the source for the data. 
The latitude and longitude were given to a tenth of a minute. 

The data on magnetic tape from NODC, EPA and Texas 
Natural Resources, required special processing due to limitations 
on processing resources and the wealth of data measurements 
recorded over relatively small areas. The NODC tape data was 
punched one card per measurement. The EPA "Storet" tape gave 
monthly data covering a span of three years. These data were 
averaged according to monthly means per station. The Texas 
Natural Resources tape was averaged by predetermined optical 
seasons (detailed in Chapter IV.C) per station. 

These data cards were then used as input for the plotting 
program. A copy of this program is included in Appendix B. 
The data points were plotted on a Mercator projection at the 
Same scale as available charts of the area. The plots were 
then overlayed onto the charts to trace the coastlines and 
bathymetric contours. Because of the irregular spacing 
of the data, a spatial averaging routine was used for areas 
where a dense Spacing of points would have caused values to 
Overprint one another. This was accomplished by a straight 
averaging of all points which were spaced within one plotted 
number width of one another. Contours were developed for 
plots by using the smallest plottable number, which was .04 


inches for the NPS Versatec Plotter. These contours were 
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then transferred using a light table to a same scale plot which 
used a number size of .10 inches, the minimum reproducible size 
for thesis presentation. Therefore, the data contours represent 
an averaging of values plotted within 1.9 miles of one another 
O he West Coast, and within 3.4 miles on the Gulf Coast. 
However, data values on the final plots represent an averaging 
of values plotted within 4.7 miles of one another on the West 
Coast and within 8.5 miles on the Gulf Coast plots. No averaging 
was performed for individual data points positioned beyond 
these limits. Figures 14 through 16 show the locations of all 
the data points used in the plots. 

For some months data were so sparse that meaningful plots 
could be obtained only by dividing the year into seasons. 
The determination of these seasons is the topic of the next 


chapter. 
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IV. TEMPORAL VARIATION OF OPTICAL PROPERTIES BY REGION 


A. INTRODUCTION 

The months of lead time required for laser bathymetric 
mission planning make it important to determine which times 
of the year provide optimum water clarity. The sparcity of 
data for certain months leads to a grouping of months into 
seasons. 

Seasonal trends were estimated by comparing the month-to- 
month variations of oceanographic factors known to affect 
turbidity with the optical measurements of those few stations 
which provided monthly optical data. Oceanographic parameters 
which correlated best with the optical measurements were used 
to delimit the seasonal boundaries during those months with in- 
sufficient data. Secular variations were not studied due to 
the lack of repetitive optical measurements at the same location 
Over a number of years. 

Coastal water turbidity can be influenced by physical, 
Chemical, and biological processes occurring both in the water 
column and in adjacent land areas. A major contributor to 
turbidity is particulate matter produced by land runoff and 
plankton, especially in areas of upwelling. Phytoplankton 
blooms can produce sudden increases in turbidity, which may be 
Closely related to upwelling which in turn may be seasonal. 

Sea surface temperature measurements can be useful in identify- 
ing areas of upwelling, and salinity measurements have been 


used to trace the outflow from large rivers such as the Columbia 
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River (Pruter and Alverson 1972). York in (1974) demonstrated 

an inverse relationship between Secchi depths and Bakun's (1975) 
upwelling index. This upwelling index was helpful in determining 
optical seasons for the west coast the U.S. Other factors 
such as large storms may increase turbidity over shorter time 
frames. 

Data for the northwest coastal waters of the U.S. were 
assembled for the region from 42°N to 49°N out to a depth of 
200 m. The California coast was not covered due to the narrowness 
of the Continental Shelf there.* 

Data from the Gulf of Mexico were divided into two groups, 
one for the area west of 89° 30'W and the other to the east of 
that line. The Western Gulf region comprised that area along 
E... of Texas and Louisiana out to a depth of 200 m. The 
Eastern Gulf region started at 24°N, 80°W and included the 
southern and western coastsof Florida and the coasts of Alabama 


and Mississippi out to a depth of 200 m. 


B. NORTHWEST COAST 

Discharge from the Columbia and other rivers and seasonal 
upwelling along the coast dominate the turbidity observed along 
the northwest coast of the U.S. Upwelling affects most of the 
coast during the summer and the river discharges produce turbidity 
plumes which affect the coastal areas near river mouths. Brown 
(1973) correlated Secchi depths for the northwest coast with 


Other simultaneously measured oceanographic parameters and found 


*For the California coast there are a number of optical obser- 
vations to which the reader can refer. References to California 
Mata which have either beam attenuation or diffuse attenuation 
measurements may be found in reports by the Allan Hancock Foundation 
(1963); Drake (1972); Frederick (1970); Karl (1976); and Winzler 
and Kelly (1977). 
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coefficients of .610 correlation with water depth, -.598 with 
Miserved Forel-Ule color codes, -.591 with Si04 and .504 with 
Salinity. Other parameters such as surface temperature, density, 
and oxygen levels did not show as good correlations. 

Bakun (1975) has attempted to quantify upwelling through 
the use of an upwelling index which is an estimate of the com- 
ponent of computed Ekman transport directed offshore. This 
transport is calculated from daily mean surface atmospheric 
pressure data. Figure 13 (Bakun, 1975) shows the upwelling 
indices computed by month versus latitude for the western 
Coast of the U.S. Shaded areas represent upwelling in units of 
cubic meters per second per 100 m of coastline. Negative values 
indicate downwelling. The figure was compiled from an average 
of 20 years of wind observations. Upwelling occurs as a result 
of change in wind patterns from southwesterly in winter to 
northerly in summer. Although the upwelling is not constant 
but varies with wind variations, this same annual cycle of 
Summer upwelling has been observed over the course of many 
years. 

The Columbia River is the only major river source of 
turbidity along the Northwest coast. It ranks second among 
Meo. rivers in volume cf discharge at its mouth with an average 
of 640,000 cubic feet per second (or approximately 7,300 m/sec) 
(Boone 1978). Peak discharges occur from May to July with 
lowest flows from August to October. Surface currents along 
this coast flow southerly from May through September, and northerly 
from mid-November through February (U.S. Navy Hydrographic 


Office 1967). This movement pushes the plume of the Columbia 
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¡River to the south during the months of peak discharge. Most 

of the remaining streams of the Northwest coastline originate 

in the Coast Range and produce highest stream flows in February 
through March and lowest flows from August to September (Geraghty 
2573). 

The preceding oceanographic factors along with K data from 
the area (Small, 1979) are summarized in Table V. The table 
shows upwelling from May to September, also the period of peak 
discharge of the Columbia River.  Upwelling and river discharge 
act together to assure higher K values at this time. A season 
Of relatively lower K was established from October to January 
when oceanic and coastal waters meet near shore and little or 
no upwelling occurs. February through April shows variability 
with upwelling occasionally establishing itself for short 
periods. Seasonal plots, figures 17 through 36, were generated 
using this division scheme. 

Figures 17 through 25 show Secchi values obtained for this 
area. Of particular interest is the delineation of the Columbia 
River plume. Figure 21 shows the 5 m contour being extended 
south across the 46° parallel towards deeper waters. Figure 22 
shows the same contour bulge as much smaller and directed northerly 
with the current at that time of year. Figure 20 shows the 5m 
our bulging straight out from the river mouth during the 
time of year when the surface currents are weak and variable in 
direction. 

Figure 26 shows all the beam attenuation measurements obtained 
for the Northwest coast for all seasons. Figure 27 shows all 


ШИЕ irradiance measurements obtained for all seasons. 
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Figures 28 through 36 show all measurements after their 
conversion to irradiance K according to the methods detailed 
in Chapter 3. Due to the abundance of Secchi measurements as 
opposed to any other type, these plots mainly reflect the 


trends indicated in the raw Secchi plots. 


ШОО WESTERN GULF OF MEXICO 

Winter storms and their associated winds and resulting 
high seas appear to be the most efficient mechanism in 
generating turbidity in the Western Gulf. Besides being able 
to resuspend bottom sediment, these storms apparently cause 
efficient tidal flushing of local estuaries and lagoons as 
suggested by Shideler (1979). These lagoons are laden with 
sediments from local rivers. However, peak discharges of 
turbid water from these lagoons is not related to high stream 
flows. Instead, the majority of sediment is trapped by the 
lagoons until storm conditions can aid high tides in dis- 
EDSrging it. 

The discharges of the Mississippi and Atchaflaya Rivers 
dominate the turbidity regime along the Louisiana Coast. 
Peak discharges for both rivers are from March to May, with 
lowest stream flows from September to November (Perret, et al 
1971). These flows dominate the nearshore coast of Louisiana 
with diminishing effect towards the outer continental shelf. 
Other stream flows show peak discharges from January to May 
and lowest stream flows from August to October. 

Coastal currents are generally weak (less than 1 m/sec 


Ever 95% of the time) but do play a role in directing turbidity 
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flows. From October to June the longshore currents are from 
the east and flow west, then south, along the Texas Coast (U.S. 
Naval Oceanographic Office, 1972). This causes the Atchaflaya 
and most of the Mississippi water to follow the Louisiana coast 
and keep turbidities high in this area. These currents also 
direct discharges from Texas lagoons to follow nearshore 
shallow waters. From July to September a longshore current 
from the south invades the southern Texas coast until it meets 
the previously described current at about 29° latitude. In- 
cursions of the Loop Current of the Gulf of Mexico onto the 
Continental shelf have been observed (Shideler, 1979). This 
shoreward incursion of deep water confines turbid coastal 
waters to the shallow inner shelf and results in lower turbidities 
elsewhere. 

Upwelling, which plays an important role in turbidity levels 
for the Northwest coast, is not thought to be important for the 
Western Gulf. Winds from the south to southwest (for Texas coast) 
and from the west (for the Louisiana coast) which would be 
needed to produce an offshore Ekman transport are infrequent 
and generally weak. 

Berryhill, et al. (1976) report finding a prevalent two-layer 
turbidity structure in their study of the South Texas outer 
continental shelf. This structure consisted of a nepheloid layer 
below a less turbid layer. This nepheloid layer varied in both 
thickness and distribution but in general became thickest toward 
the outer shelf. This type of structure casts doubt on the 
validity of the use of Secchi measurements to obtain K for the 
entire water column because they only indicate transparency 


near the surface. 
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Table VI illustrates the variables used in determining 
optical seasons. Wind data were obtained from a Summary of 
Synoptic Meteorological Observations (U.S. Naval Weather 
Service Command, 1970). Sea heights were obtained from a 
U.S. Naval Oceanographic Atlas (1972). Zooplankton counts 
were obtained from a study done by the Louisiana Wild Life 
and Fisheries Commission (Perret, 1971). Secchi transects 
were from a study of the South Texas outer continental shelf 
(Kamykowski, 1979). 

A season from November to February was chosen to represent 
that time of year when storms are most prevalent and an efficient 
tidal flushing of lagoons is thought to occur. There is a varia- 
ble season from March to May when winds are occasionally strong 
and outflows from the larger streams are high. June to October 
is the third season, when winds and streamflow are lowest. This 
third season is thought to consistently provide the best overall 
water visibility. Hurricanes which occur during this season 
could raise turbidities substantially, but are infrequent and 
usually short term in their effects. 

Figures 37 through 39 show Secchi measurements made in the 
Western Gulf. Within the lagoons of Texas and the nearshore 
area of Louisiana, visibility is most related to streamflow 
and therefore March to May is generally the period of shallowest 
Secchi depths and when streamflows are highest. June to October 
is generally that period when the deepest Secchi depths were 
found, but a sparcity of data during these months makes compari- 
ΠΝ. citficult. November through February generally are the 
months when highest turbidities are encountered in the offshore 


deeper waters. 
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Figures 40 and 41 show the only beam attenuation data col- 
lected for this area. The more turbid waters of November to 
February of the two figures provide a further indication of 
tidal flushing of lagoons during this stormy season. 

Figures 43 through 45 show the conversion of all the data 
to the common parameter K. Secchi depths less than 2m werte not 
dealt with by the 4 to K conversion formula which was derived 
from measurements deeper than 2 m. Therefore, the Texas 
lagoon measurements do not appear on these figures. These 
plots display the same seasonal variations mentioned before. 
Shideler (1975) notes that secular variations were thought to 
be related mainly to annual stream flow variations and storm 


occurrences for those years. 


ШІ EASTERN GULF OF MEXICO 

The eastern Gulf data indicate similar mechanisms of 
turbidity generation as were present in the Western Gulf. The 
main differences are due to the appearance of some upwelling 
during winter months and a change of the months of highest 
stream flow along the Florida coast. 

Alexander, et al. (1977) attribute mixing of the water 
columnover theeasternGulf'scontinental shelf to wind stress 
which can be linked to the frequency of low pressure disturbances 
(storms) crossing this area. In winter when storms are most 
frequent the water column was observed to be unstable and well 
mixed. By contrast, in summer and fall the water column is 
Stable with established thermoclines and haloclines. 


Alexander, et al. (1977) also noted that the passage of Hurricane 
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Bloise did not raise turbidity levels to the degree that per- 
sistent winter storms had. Also, the rise in turbidity caused 
née hurricane disappeared after one week. 

Monthly Ekman transports for 2 degree squares in the 
Eastern Gulf have been computed by Ichiye, et al. (1973). Their 
results show transports to be lowest from June to August with 
averages less than 3,000 gm/cm-sec. Transport is generally 
alongshore except from December to February when offshore 
transports result from scattered areas of upwelling along the 
Florida coast, especially near De Soto Canyon. In other months 
there may be occasional periods of offshore transport but of 
smaller magnitudes. 

Surface currents generally flow in a northwesterly direction 
along the coast in the nearshore areas. This directs Mississippi 
River flow away from the Eastern Gulf except from April to 
September when a reversal in flow direction occurs from the 
Mississippi River mouth to the Florida panhandle.  Incursions 
of pockets of Mississippi outflow were observed by Alexander, 
et al. (1977) along with Loop Current eddies which had intruded 
onto the continental shelf during these summer months. No Loop 
Current intrusions were found during other seasons. 

Table VII summarizes the main parameters used to delineate 
optical seasons. Wind mixing Occurring during winter storms 
and a simultaneous increase in upwelling makes January to March 
a season of highest K. Lowest K values can be expected from 
June to September when the water column is most stable and 
clear water from the Loop Current intrudes upon the outer 


shelf. Exceptions to this occur along the Mississippi, Alabama, 
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and western Florida coasts where turbid Mississippi River out- 

flows make K higher. The two remaining seasons are transition 

periods where K is variable and less predictable. Stream flows 
other than the Mississippi do not appear to influence turbidity 
levels greatly. 

Eres 46 through 53 show the Secchi measurements obtained 
for the Eastern Gulf. Figures 54 through 59 show beam and 
irradiance attenuation measurements and figures 60 through 66 
show all measurements converted to irradiance K. Overall, 


these plots show the lowest K values encountered in this study. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


А. CONCLUSIONS 

The present and near-future water penetration capability 
of laser systems is limited to about four attenuation lengths 
(Enabnit, 1979) as was discussed in Chapter II (see also Fig. 6). 
Seasonal Kd values were plotted in nineteen figures for Pacific 
Northwest and Gulf coastal areas. These figures indicate where 
and when laser bathymetry systems would be successful, given 
ENSEROS4 criterion. These seasonal Kd values were derived from 
Optical data collected over different years. However, inter- 
annual variations may change the utility of laser bathymetry 
systems in a particular area. 

In the coastal areas of Oregon and Washington laser bathy- 
metry would probably not be successful at any time over signifi- 
cant areas as is shown in Figures 67 to 75. A Kd value of ten 
was contoured (the lowest contourable value) as an indicator 
of possible areas surveyable by laser. The area within the 
Kd-10 contour is small for all seasons; the even smaller area 
Ж ЕЕЕ Kds4-is too small for practical laser hydrographic opera- 
tions. 

In the Western Gulf, limited areas are candidates for laser 
bathymetry as shown by Figures 76 to 78. The southern half of 
Laguna Madre, Texas, is survevable by laser bathymetry from 
November through February. The area from Matagorda Bay to 
Sabine Pass off Texas is surveyable all year up to depths of 


10 to 20 fathoms as indicated by the Kd=4 contours. West of the 
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Mississippi Delta, possibly to Sabine Pass, a strip from outside 
local estuaries and bays to the 10-fm depth contour is available 
ШО laser bathymetry from March to October. 

The Eastern Gulf coastalarea is the area best suited to 
laser bathymetry as indicated in Figures 79 to 85 by the Kd=4 
contour. Off Florida from Panama City to the Florida Keys, an 
area of 30,000 nmi? bordered by the 30-fm depth contour may be 
surveyed by laser bathymetry during June through September, 
and a reduced area bordered by the 20-fm depth contour all year. 
This 30,000 nmi* area alone represents 15 years of work for one 
laser system and represents a costs savings of $69 million (1977 
dollars) over launch/sonar surveys. From Panama City to the 
eastern Mississippi Delta an area of 8,800 плі“ bordered by 


the 20-fm depth contour is surveyable from October through 


December. 


B. RECOMMENDATIONS 

l. More optical measurements, especially irradiance and 
beam attenuation data, should be collected throughout 
the year in U.S. coastal waters. 

2. More simultaneous measurements of Z! K, and C at the 
same wavelength should be taken and more investigations 
of the relationships between them should be made. 
Universal relationships should not be expected due to 
regional variation in suspended particies and yellow 
substance (Gordon and Wouters, 1978) and, therefore, 


such studies should be regional. 
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Investigators should send optical measurements to the 
National Oceanographic Data Center, especially measure- 
ments made in the Gulf Coast area. NODC should establish 
irradiance and beam attenuation data categories as part 
of the Ocean Stations Data Base. 

This thesis studied seasonal variation of marine optical 
properties in U.S. coastal waters; others should investi- 
gate long term trends and their cause. 

Investigators should always state the wavelength at which 
their optical data were taken. 

Prior to the commencement of laser hydrographic survey 
operations, Nimbus 7-CZCS data for the areas of interest 
should be examined and used to update the figures pre- 
sented here. 

This thesis did not examine the bottom reflectivities of 
coastal waters because of their generally relatively minor 
effect compared to attenuation lengths as a limiting 
parameter for laser bathymetry operations. However, in 
areas of marginal utility for laser bathymetry, bottom 
reflectivity may become significant. This suggests 
further study of this effect is needed. 

The National Ocean Survey should consider the west coast 
of Florida for airborne laser hydrographic operations 


due to the favorable marine optical environment there. 
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APPENDIX A 


DENE SOUREERZOFZTN-STTÜU OPTICAL MEASUREMENTS 
IN COASTAL WATERS 

Data Source: CDS Data Base 
Institution: Texas Natural Resources Information System 
Geographic Area: Texas 
Date of Tape Generation: May 1979 
Data Range: 1968-77 
Optical Data Type: 2 


5 
Number of Data Observations: 14,232 
Шаға Reduction Procedures: Data was on magnetic tape. 


Converted to K by methods of Section III.B.3. 
Data Source: SD-2 Oceanographic Stations (Master Records) 
Data Base 
Institution: National Oceanographic Data Center 
Geographic Area: Primarily Oregon and Washington 
Date of Tape Generation: March 1979 
Data Range: 1952-74 
Optical Data Type: 2- 
Number of Data Observations: 1329 
Data Reduction Procedures: Data was on magnetic tape. 
Converted to K by methods of Section III.B.3. 
Data Source: STORET Data Base 
Institution: Environmental Protection Agency 
Geographic Area: Florida 
Date of Tape Generation: June 1979 
Mata Range: 1966-79 
Optical Data Type: 2- 
Number of Data Observations: 3,000 
Data Reduction Procedures: Data was on magnetic tape for 
selected stations off Florida. Converted to K by methods 
ersSection III.B.3. 
Investigator: Barret, Barney B. 
Data Source: Barret (1971) 
Institution: Louisiana Department of Wildlife and Fisheries 


Geographic Area: Louisiana 
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Date of Data Collection: 1968-1969 


Optical Data Type: Z 


Data Reduction Procedure: Computed K from Z 
emteiined in Section III.B.3. 


= ру method 
Investigators: Barret, Barney B. et al 

Data Source: Barret et al (1978) 
Institution: Louisiana Department of Wildlife and Fisheries 
Geographic Area: Louisiana 

Date of Data Collection: 1974-76 

Optical Data Type: 2 


5 


Data Reduction Procedure: Computed K from 2 。 by method outlined 
9 Section Ill1.B.3. 

Investigators: Berryhill, Henry L., Jr., et al 

Data Sources: Berryhill et al (1976) and Berryhill et al (1978) 

۶۰۰٣۰٥0105: U.S. Geological Survey 

Geographic Area: Texas 

Date of Data Collection: 1975-1977 

Optical Data Туре: С 

Instrument: Martek Transmissometer (528nm) 

Data Reduction Procedure: Integrated C computed by method 
ЫШ л леа in Section III.B.2 for 28 stations sampled six 
different times. Most stations had the beam transmission 
measurements (surface, middle, and bottom) per profile but 
some lacked bottom measurements. For these stations the 
middle value was used to the bottom depth. K was then 
computed from C. 

Mavestigators: Carder, K. L., and Haddad, K. D. 

Data Source: Carder and Haddad (1979) 

Me titution: University of South Florida 

Geographic Area: Mississippi, Alabama, and Florida 


Date of Data Collection: 1976-1978 


Б, 





Optical Data Type: C 

Instrument: Hydro Products Transmissometer (550nm) 

Data Reduction Procedure: The original data were given as 
suspended particles beam attenuation coefficient (C_) 
contours on profile figures. The use of this data 
required special processing. First, values of C_ had to be 
picked off the figures for the C_ vs depth profile at each 
ocean station. Second, C_ was cónverted to C for each layer 
by using the relation C_ + .06935 = C (Carder and Haddad, 
eo). Third, the methBds of Section IIl.B.2. were used to 
obtain K. 

Investigator: El-Sayed, Sayed 27. 

Data Source: Unpublished Data (El-Sayed 1974) 

Institution: Texas A&M University 

Geographic Area: Gulf of Mexico 

Date of Data Collection: 1971-73 

Optical Data Type: 2 


Data Reduction Procedure: Computed K from 2. by method 
ИШ леа in Section III.B.3. 

Investigators: Godcharles, Mark F. and Jaap, Walter C. 
Data Source: Godcharles and Jaap (1973) 


Institution: Florida Department of Natural Resources Marine 
Research Laboratory 


Geographic Area: Florida 
Date of Data Collection: 1970-71 
Optical Data Type: 2- 


Data Reduction Procedure: Computed K from 2- by method 
u lined in Section III.B.3. 

Investigators: Gordon, Howard R. and Dera Jerzy 

Data Source: Gordon and Dera (1969) 

Institution: Institute of Marine Sciences, University of 


Miami 


52 





Geographic Area: Florida 
Date of Data Collection: December 1967 
ШШ са! Data Type: K 


Instrument:  Irradiance meter with a selenium photocell 
covered by a 525nm filter and a diffuse screen. 


Data Reduction Procedure: None 


Investigator: Grady, John R. 

Data Source: Grady (1979). Five Cruises of GUS III 
Institution: National Marine Fisheries Service 
Geographic Area: Louisiana and Texas 

Date of Data Collection: Jan-May 1966 

Optical Data Type: 2 


S 


Data Reduction Procedures: Z converted to K by methods 
a Section III.B.3. 

Investigators: James, W. P., et al. 

Data Source: James (1977) 

Institution: Texas A&M University 

Geographic Area: Florida 

Ele Of Data Collection: 2/76, 3/76, 7/76 

Metical Data Туре: К, 2 


5 


Instrument: Kahl Scientific Co. Universal Radiometric 
Submarine Photometer (Green filter) 


Data Reduction Procedures: К data reduced to a depth_ averaged 


method of Section II1.B.l. Z data converted to K by 
method of Section III.B.3. 


Investigatars: Joyce, Edwin A. and Williams, Jean 
Data Source: Joyce and Williams (1969) 
Institution: Florida Department of Natural Resources 


Geographic Area: Florida 
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Date of Data Collection: 1966-69 
Optical Data Type: 2 


Nasa Reduction Procedure: Ὃ عون‎ K from 2. by method 
outlined in Section III.B. 

Investigators: Kamykowski, Daniel et al. 

Data Source: Kamykowski et al (1978) 


Institutions: University of Texas Marine Science Institute 
Port Aransas Marine Laboratory 


Geographic Area: South Texas 
Date of Data Collection: 6/718-168 
Optical Data Type: Quantum K, C, and Ze 


Instrument: Lambda Photometer (quantum) 
Martek Transmissometer (528nm) 


Ша Reduction Procedure: Quantum K reduced to K by procedure 
outlined in Section 111.B.1. Integrated C computed by method 
outlined in Section III.B.2. nen K_ and Z_ were used to 
compute the conversion formula, to К for the Western Gulf 


Coast and the procedure is E e 1707-560617 01 ۲٢76:8 3 All 
data for one oceanographic station. 


Investigators: Manheim, Frank T.; Steward, Robert G., and 
Carder, Kendall L. 

Data Source: Manheim, Steward and Carder (1977). 
Institution: University of South Florida 

Geographic Area: Mississippi, Alabama, and Florida 

шіге OL Data Collection: 1975-1976 

Optical Data Type: C 

Instrument: Hydro Products Transmissometer (550nm) 

Data Reduction Procedures: C and K computed by methods of 
Section III.B.2. 

Investigators: McGrail, David W.; Huff, David; Jenkins, Stacy 
Data Source: McGrail, Huff and Jenkins (1978) 

Institution: Texas A&M University 
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I. 


۰. 


Geographic Area: Texas and Louisiana 

Date of Data Collection: 1977 

Optical Data Type: C 

Instruments: Martek Transmissometer 

Data Reduction Procedures: The data published were trans- 


mittance profiles graphed for each ocean station. Thus, 
inflection points were picked off and the methods of Section 


III.B.2 were used to computer C and К. 


Investigator: Oregon State University Cruises 
Data Source: Unpublished data (Pak 1974) 
Institution: Oregon State University 
Geographic Area: Oregon 

meee Of Data Collection: 1968-1971. 

Optical Data Type: 27 


5 


Date Reduction Procedure: Computed K from Z by method 
шише ппей in Section III.B.3. 
Investigators: Pak, Hasong and Zaneveld, Ronald V. 


Date Source: Unpublished data (Pak 1979) Conclusions of 
this data were published by Pak and Zaneveld (1977). 


Institution: Oregon State University 

Geographic Area: Oregon 

Date of Data Collection: 8/74-5/75 

Optical Data Type: C 

Instrument: OSU Transmissometer (660nm) 

Data Reduction Procedure: Computed C at 660nm converted to 
C (532nm), and converted to K by methods outlined in Section 
I1 85.2. 

Investigator:  Saloman, Carl H. 

Data Source: Saloman (1974) 


Institution: National Marine Fisheries Service 


Geographic Area: Florida 
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Date of Data Collection: 11/70-9/71 
Optical Data Type: Z 


Data Reduction Procedure: Computed K from Z_ by method out- 
mea in Section III.B.3. 7 

Investigators: Saloman, Carl H. and Collins, L. Alan 

Data Source: Saloman and Collins (1974) 

Institution: National Marine Fisheries Service 

Geographic Area: Tampa Bay, Florida 

Date of Data Collection: 1/71-12/72 

Optical Data Type: 2 


S 


Date Reduction Procedure:  Computed K from 2. by method 

Buelined in Section III.B.3. 

Investigators: Small, Lawrence F. and Curl, Herbert, Jr. 

Data Source: Unpublished Data (Small 1979). Conclusions 
of this data were published by Small and Curl 
(1968). 

Institution: Oregon State University 

Geographic Area: Oregon 

Date of Data Collection: 4/62-4/65 

Optical Data Type: K 


Instrument: Kahl Scientific Irradiance Meter (photopic) 


Data Reduction Procedure: Mne 


Investigators: Stevenson, W. H., and Pastula, E. J. 
Data Source: Stevenson and Pastula (1973) 
Institution: National Marine Fisheries Service 
Geographic Area: Mississippi Sound 

Date of Data Collection: August 1972 

Optical Data Type: 2 


S 


Data Reduction Procedure: Computed K from Z by method out- 
fined in Section III.B.3. 
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TABLE ІІ. DATA USED TO CONVERT SECCHI DEPTH TO IRRADIANCE 
ATTENUATION (Z TO к.) FOR EASTERN PACIFIC COASTAL 


WATERS 
1957 North North 5 
Station Date Latitude Longitude Z. (m) Ko (m^!) 
49A 7/9 39° 59° 1269 38' 32 .0718 
5] 7/10 409 02' 127° 07 32 .0363 
61 7/15 44° 50' ος, ΙΙ; 22 1114 
64 7/23 πο ης) 126% 38' 20 A 
65 7/23 43° 51' 02000 32 .0708 
66 7/24 430 43: П 32 - .0747 
67 ДЕ 15013 πο 26 0903 
83 7/30 390 01' 1267 26' 28 .0674 
88 8/4 38% 56' 125% 42' 26 .0651 
89 8/4 389 57' 126% 10' 24 .0900 
9] 8/5 38% 14' 1250 40! 20 .1096 
92 85 38° 05' 1259 25" 16 .1327 
93 8/5 Ba 125° 44: 22 ‚0938 
95 8/6 37% 06' 125% 07' 24 ‚0999 
96 8/6 369 58' 1247 56' 18 1291 
98 8/6 369 45" 1242 ται 20 .1340 
101 8/7 36° 45' едал“ 22 m 
102 8/8 36° 08' lOS e 20 „1151 
103 8/8 35% 19' ο 10 ‚2577 


Original station data from Callaway and McGary (1959); 2 computed 
by van Norden. 
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TABLE III. 


Date 


6/29/78 
6/29/78 
7/24/78 
7/25/78 
7/25/78 
9/25/78 
9/26/78 
11/08/78 
11/08/78 
11/09/78 
11/09/78 


All observations at N 229 34' Latitude and W 96° 50! 


Longitude. 


obtained from Jerlov 


Local 
Time 


1300 
1700 
1200 
0800 
1200 
1345 
0000 
1200 
1600 
0800 
1200 


Z, (m) 


23815 


20 
26 
19 
20 


14 


Quanta К. 
à э 


.0609 
20903 
.0683 
. 7 
.0442 
ο 
1295 
5089 
192 
.187 
.158 


K 7 


Q 


1.085 


1.06 
1.09 
1298 
1212 
ΙΟ 
975 


E 


ου 


.980 
1990 
5950 


K 


s m 


.0661 


.0964 


‚0748 


«0634 


‚0495 


ην 


209 


ди 


‚180 


.177 


.143 


DATZEUSEDZIOZEONVERT SECCHI DEPTH TO IRRADIANCE 
ATTENUATION (2, TO К.) FOR WESTERN GULF COASTAL 
WATERS 


pis 


Original station data from Kamykowski, et al. 
89798); Quanta Ko and K computed by van Norden and K/ Kg 


(05976) 
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FIGURE 1. LASER/SONAR COST COMPARISON, COSTS IN 1977 DOLLARS 
ШЕС, 1979)” 
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SPECTRAL IRRADIANCE ATTENUATION COEFFICIENTS OF 
DOWNWARD IRRADIANCE IN THE SURFACE LAYER FOR 
JERLOV OPTICAL WATER TYPES (AVCO EVERETT RESEARCH 
LABORATORY, INC., 1978). 
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` FIGURE 5. SYSTEM ATTENUATION COEFFICIENTS AS A 


FUNCTION OF RECEIVER FOV AND AIRCRAFT 
ALTITUDE. DATA FROM KRUMBOLTZ (1979). 
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ATTENUATION LENGTHS TO THE BOTTOM. 
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PERCENTAGE OF SURFACE QUANTA (350-700nm) 





FIGURE 7 


ШИ ROriILES OF PERCENTAGE OF SURFACE 
QUANTA(350-700nm) FOR DIFFERENT WATER TYPES 
ειν", 7ο, 
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FIGURE 8. THE RATIO OF IRRADIANCE TO QUANTUM IRRADIANCE 
IN THE SPECTRAL RANGE 350-700 nm AS A FUNCTION 
OF DEPTH IN DIFFERENT OPTICAL WATER MASSES 
MOS 76) 


15 





K,-1.2525'«0.05 


22-078К--004 


IRRADIANCE ATTENUATION COEFFICIENT K (m7) 





0 1 ‚2 3 4 5 
SECCHI DEPTH RECIPROCAL 22(М71) 


FIGURE 9. LINEAR REGRESSIONS OF E AND Z. FOR EASTERN PACIFIC 
COASTAL WATERS. 
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FIGURE 10. LINEAR REGRESSIONS OF κ. AND ὃς FOR WESTERN 
GULF COASTAL WATERS. 
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FIGURE 11. LINEAR REGRESSIONS OF K AND 2 FOR POOLE AND 
ATKINS DATA. 
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WITT,ET AL(1976) 
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FIGURE 12. LINEAR REGRESSIONS USED TO CONVERT 2. TO 
p FOR SELECTED U.S. COASTAL WATERS. 
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